First-principles investigation of spin polarized conductance in atomic carbon wires 
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We analyze spin-dependent energetics and conductance for one dimensional (ID) atomic carbon 
wires consisting of terminal magnetic (Co) and interior nonmagnetic (C) atoms sandwiched between 
gold electrodes, obtained employing first-principles gradient corrected density functional theory and 
Landauer's formalism for conductance. Wires containing an even number of interior carbon atoms 
are found to be acetylenic with a — n bonding patterns, while cumulene structures are seen in wires 
containing odd number of interior carbon atoms, as a result of strong vr-conjugation. Ground states 
of carbon wires containing up to 13 C atoms are found to have anti-parallel spin configurations of the 
two terminal Co atoms, while the 14 C wire has a parallel Co spin configuration in the ground state. 
The stability of the anti-ferromagnetic state in the wires is ascribed to a super-exchange effect. For 
the cumulenic wires this effect is constant for all wire lengths. For the acetylenic wires, the super- 
exchange effect diminishes as the wire length increases, going to zero for the atomic wire containing 
14 carbon atoms. Conductance calculations at the zero bias limit show spin-valve behavior, with the 
parallel Co spin configuration state giving higher conductance than the corresponding anti-parallel 
state, and a non-monotonic variation of conductance with the length of the wires for both spin 
configurations. 



I. INTRODUCTION 



Miniaturization from sub-micron conventional solid 
state devices to extreme small scale single organic 
molecule based devices has been the focus of inten- 
sive research in recent years, motivated by the emerg- 
ing field of molecular scale electronics and quantum in- 
formation technology. Controlled transport of electrons 
in molecular wires containing only a few atoms forms 
the basis of molecular scale electronics. Significant re- 
cent advances in experimental techniques have made it 
possible to fabricate nano-wires containing only a few 
atoms and measure their electrical properties.^ Specif- 
ically, atomic carbon wires containing up to 20 atoms 
have been synthesized.^ These carbon wires serve as ideal 
models to develop understanding of and to eventually 
control the mechanism of electron transport in finite 
one-dimensional (ID) systems. Previous theoretical and 
experimental studies'^"^ on atomic and molecular wires 
have been primarily limited to charge transport, with 
some recent exceptions. Recent experimental mea- 
surements have shown that electron spin polarization can 
persist considerably longer than charge polarization. -'^^ It 
is consequently highly desirable to learn how to manipu- 
late and enhance the control of electron transport offered 
by spin degrees of freedom, adding another dimension to 
the emerging field of molecular scale electronics and re- 
vealing important information for potential applications 



in spin-based molecular electronics (spintronics) as well 
as in quantum information processing. 

Pure carbon clusters have a long history. ^^'^^ Clusters 
with less than 10 atoms are known to have low-energy 
linear structures characterized by cumulenic bonding 
(C=C==C=C) with near equal bond lengths. These struc- 
tures are stabilized by strong 7r-conjugation between 
the double bonds which are alternately directed in the 
X- and y-planes perpendicular to the bonds. Some of 
these clusters also possess cyclic forms, which also be- 
come the stable form for larger sizes. Lang and Avouris 
calculated the conductance of such a cumulenic carbon 
atom wire, i.e., with all C-C bond lengths constrained 
to be equal, connected on both ends directly to metal 
electrodes.^ They found an oscillatory behavior, with 
wires composed of odd numbers of carbon atoms having 
a higher conductance than even-numbered wires. This 
was contrary to expectations based on a simple molec- 
ular orbital theory of the cumulene structure, and was 
attributed to electron donation from the metal contacts 
into additional 7r-bonds formed between the terminal car- 
bon atoms and the electrodes. 

The presence of terminal magnetic atoms has been re- 
cently been shown^** to modify both the structures and 
conductance properties of these wires. Pati et al have 
reported first-principles calculations^^ of spin-dependent 
electronic structures and energetics, as well as spin- 
polarized conductance of small carbon wires containing 
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up to five carbon atoms that arc terminated by mag- 
netic atoms which are in term attached to gold elec- 
trodes. The magnetic atoms can act as spin polarizers 
or filters, resulting in a strong spin valve efi'ect. These 
results showed that when terminated by magnetic atoms, 
the TT-conjugatcd structure is not necessarily the lowest 
energy structure. In particular, for the even-numbered 
carbon wires, the acetylenic structure with alternating a- 
and multiple 7r-bonds becomes more stable. The calcula- 
tions also showed that the ground states of these magnet- 
ically terminated carbon wires have anti-parallel terminal 
atom spin configurations, which could be rationalized by 
contribution from super-exchange effects. The conduc- 
tance of the wires was seen to non-monotonous with wire 
length just as in the pure carbon wires, ^ and is shown 
to depend strongly on the magnetic configuration of the 
terminal atoms, with the result that the wires could be 
made to act as a molecular spin valve. 

This first study of spin-dependent properties of short 
carbon wircs"'^^ raised a number of fundamental questions. 
In particular, what happens when the length of the car- 
bon wire spacer between the magnetic atom species is in- 
creased? How do the wire structures and the relative en- 
ergetics of the different spin configurations change? how 
do the ground state spin configurations depend upon the 
length? How does the effective exchange coupling be- 
tween the magnetic species change with the number of 
carbon atoms in the wire? What is the critical length 
of the connecting wire up to which the super-exchange 
effects identified for the short wires survive? 

In order to address these questions and to thereby 
improve our understanding of spin-dependent electron 
transfer in extended molecular systems, we extend here 
the work in ref.l4 up to wires containing 14 carbon 
atoms. Our first principles calculations, which explic- 
itly include spin-polarization effects, reveal that the anti 
parallel Co spin state remains the ground state for wires 
with up to 13 carbon atoms, and also show a continua- 
tion of the alternation between cumulcnic and acetylenic 
structures for odd and even wire lengths respectively. 
However, the 14 carbon atom wire is seen to have a 
parallel Co spin configuration in the groimd state. In- 
terestingly, we find that for wires containing an odd 
number of carbon atoms the energy difference (AE) 
between the anti-ferromagnetic and ferromagnetic state 
(anti-ferromagnetic being the ground state) remains con- 
stant as a function of wire length. In contrast, in the 
acetylenic carbon chains, this energy difference decreases 
exponentially as a function of the number of atoms with 
the exception of the 2-carbon atom wire. Analysis of 
this change in ground state spin configuration (14 carbon 
atom wire) in terms of the super-exchange contribution 
allows us to estimate the characteristic length for super- 
exchange in acetylenic carbon wires as ~20 A. We also 
find that the 7r-conjugated cumulenic wires exhibit higher 
conductance than the acetylenic wires. Finally, the cal- 
culated magneto-conductance for different wire lengths 
shows a large difference between the two magnetization 



states, particularly for C-wires containing 13 and 14 C- 
atoms, suggesting its potential applications in molecular 
magneto-electronics. 

The remainder of the paper is organized as follows. 
Our computational approach is described in Section II. 
The results and discussions are presented in Section III. 
Section IV summarizes our results. 



II. COMPUTATIONAL DETAILS 

As in the previous study of short atomic wires, -'^^ we 
utilize an architecture consisting of chains of non mag- 
netic C-atoms connecting two magnetic Co atoms. The 
Co-(C)„-Co wire structures, with n=6 to 14, are subse- 
quently inserted between two metal gold electrodes for 
calculation of spin polarized conductance. In a magnetic 
system like this, the total conductance can be evaluated 
as: 

9t ~ 9 spin— conserved ~^ 9 spin— flip j (f) 

where gspin-conserved is the Conductance from the spin 
conserved part and gspin-fHp is the conductance due to 
spin flip scattering. The latter contribution plays a sig- 
nificant role only when the spin-orbit coupling plays a 
significant role. Since the spin flip scatterthe spin-orbit 
coupling effect in highly ordered, strongly conjugated C- 
wire structures is expected to be small, leading to a large 
spin-flip scattering length, we have assumed the scatter- 
ing to be coherent and have not included relativistic spin- 
orbit coupling effects in the present paper. The spin con- 
served part of the conductance is calculated as: 

9 spin— conserved ~ 9 ~t~ 9^ •) (^) 

where g" and g'^ are the contribution to conductance 
from up (a) and down (b) spin states, respectively. Since 
at low bias the conduction primarily occurs in the close 
proximity of the Fermi energy of the metal contact, we 
can use Landauer's approach^^ to calculate and at the 
Fermi energy. In the zero bias limit, we have: 

9'^^^\Ef)= ^ T«(«(E/), (3) 

where T°'^^\E f) is the transmission fimction for the spin 
up (a) or spin down (/3) electrons. This is evaluated using 
the Green's function derived from the Kohn-Sham ma- 
trix obtained from self-consistent spin unrestricted Den- 
sity Functional calculations."'^^ We have employed a gra- 
dient corrected Perdew-Wang 91 exchange and correla- 
tion functional^^ and double numerical polarized basis 
sct^" for the calculation of energetics and magnetic struc- 
tures. Both the spin configurations and geometry for 
parallel and anti-parallel magnetic states between the 
Co atoms are simultaneously optimized using the self- 
consistent DFT approach. Anti-parallel magnetic config- 
urations between the Co atoms are obtained by making 
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use of the broken symmetry formalism. Details of this 
procedure can be found in refs 13 and 14. From the cal- 
culated spin-polarized conductance, we then estimate the 
magneto conductance (MC) according to: 



MC = 



5t(T,T) 



(4) 



where gt{'\, T) and 5t(T, J.) are given by the total conduc- 
tance, Eq.(l) in the parallel and anti-parallel configura- 
tions, respectively. 



III. RESULTS AND DISCUSSION 



A. Structures, magnetic properties, and energetics 



Using the procedures summarized in the previous sec- 
tion, we have optimized both the spin state and geom- 
etry for the Co-(C)„=6-i4-Co wire structures in parallel 
and anti parallel magnetic configurations between the Co 
atoms. Similar to the earlier report for short atomic wires 
containing up to 5 carbon atoms. we find a clear tt — tt 
and a — IT bonding pattern for the C-atoms in the wire. 
This is illustrated in Fig. 1 for n=ll and 12 carbon atom 
wires, which shows the ground state 
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structures of both magnetic configurations. The bond 
distances in the even carbon wires show a clear alterna- 
tion, both for ferro and antiferro, consistent with a — tt 
bonding, antiferro, consistent with a — tt bonding. In 
contrast, the odd wire show no evidence of bond al- 
ternation, consistent with pure 7r-conjugated structure. 
Comparing the energy for parallel(ferro) and anti par- 
allel(antiferro) spin configurations in the wires, we find 
that the anti-ferromagnetic state is lower in energy for 
all wires studied here, except n=14, which shows a ferro- 
magnetic ground state. The calculated energy difference 
between ferromagnetic and anti ferromagnetic configura- 
tions, AE = E{1, 1) — E{'\, 1), is shown as a function of 
number of C-atoms in the atomic wire in Fig. 2. For com- 
parison purposes we have also shown here the results for 
short atomic wires obtained in Ref.lA. The energy dif- 
ference between the two magnetization states are found 
to be larger than fcsT at room temperature, suggesting 
that these antiferro-magnetic states are stable in normal 
operating conditions. The lower energy for the anti fer- 
romagnetic spin configuration between the terminal Co 
atoms can be attributed to a super-exchange interaction 
that is facilitated by strong overlap of the magnetic Co 
and the non-magnetic C-atoms. A careful analysis of AE 
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FIG. 1. Schematic of the bonding pattern in a carbon wire 
containing (a) 11 and (b) 12- carbon atoms. Terminal black 
spheres are Cobalt atoms and gray spheres are carbons. All 
the optimized bond distances are indicated for both paral- 
lel(ferro) and anti-parallel (antiferro) spin configurations of 
the atomic carbon wire. 
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FIG. 2. The energy difference {AE = E{1, i) - E{t T)) 
between parallel and anti parallel spin configurations in car- 
bon atom wires, as a function of the number of carbon atoms 
in the wire, n. Even n chains are connected with a dashed 
line and odd n chains are connected with a solid line. Values 
for n=l-5 carbon atoms are taken from Ref.lA. The energy 
difference(Ai?) is given in eV. 

as function of wire length suggests that for the strongly 
TT-conjugated wires i.e. cumulenic structures (odd num- 
ber of C-atoms), the energy difference is approximately 
independent of the number of C-atoms in the wire. For 
CT — vr conjugated C-wires (even number of C-atoms), AE 
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exhibits an non-monotonic behavior with wire length. In 
particular, the anti-ferromagnetic state for C-wires con- 
taining 2, 4 and 6 carbon atoms are more stable than 
those for the 7r-conjugated C-wires. For acetylenic wire, 
AE decreases in an exponential manner (with the ex- 
ception of the 2 carbon atom wire) and is found to be 
negative for the wire containing 14 C-atoms. This sug- 
gests that the super exchange, which stabilizes the anti- 
ferromagnetic phase in cr — TT conjugated system, atten- 
uates exponentially and becomes negligible for a wire 
containing 14 C-atoms. This allows us to estimate the 
super-exchange characteristic length for a cr — tt conju- 
gated system to be ~20 A. 

The additional stability for the wires containing 2, 4 
and 6 C-atoms compared to 1, 3 and 5 C-atom wires 
could be explained as follows. In wires containing 1, 3 
and 5 C-atoms, the exchange interaction is facilitated 
between the two terminal Co atoms through delocalized 
spins shared by both the Co atoms stabilizing further the 
ferromagnetic coupling compared to that in short even 
atom wire. Similar ferromagnetic stabilization has been 
seen in the Fe^+-Fe^+ compounds. This extra stabil- 
ity of ferromagnetic ordering in the short odd atom wire 
leads to smaller AE as compared to wires containing 2, 
4 and 6 C-atoms. Thus both double exchange and super 
exchange play an important role in stabilizing the mag- 
netic ordering in these systems. In odd C-wires, dou- 
ble exchange and super exchange effects remain constant 
due to continuous p-conjugation. However, the super ex- 
change effect exceeds the double exchange effect result- 
ing in anti ferromagnetic configurations as the ground 
states. In even C-wires containing 2,4 and 6 C-atoms, 
the double exchange effect is smaller compared to 1, 3, 
and 5 C-wires, respectively. This is due to the lack of 
spin delocalization (Fig. 3) in even wires, which in effect 
destabilizes ferromagnetic ordering leading to a large en- 
ergy difference between the two magnetic configurations 
(A^ in Fig. 2). 
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FIG. 3. Atomic spin density of the terminal Co atoms for 
both parallel and anti-parallel Co spin configurations, shown 
as a function of the number of carbon atoms in the wire, n. 



B. Spin Polarized Conductance. 



Using Landauer's approach as outlined in Section II, 
we have calculated the spin-polarized conductance in the 
zero bias limit. The results are summarized in Fig. 4. 
Several interesting features are apparent here. First, the 
conductance in the C-wire is found to be higher for par- 
allel than for anti-parallel spin configuration of the ter- 
minal Co-atoms. This is a prerequisite for the spin-valve 
effect, which is primarily due to spin dependent scat- 
tering and which has been observed in magnetic/non- 
magnetic hetero bulk structures. ^^'■^■^ Second, both paral- 
lel and anti-parallel spin configurations show oscillations 
in conductance as a function of the wire length. For the 
parallel spin configurations, the conductance oscillation 
is damped after n=8 C-atoms and remains almost con- 
stant at about Igo (ffo = 2e'^/h) for the wire with n=12, 
13 and 14 C-atoms. In contrast, in the anti-parallel case, 
the conductance is seen to decrease as n increases and to 
finally vanish for a wire for 12 and 14 C-atoms. The faster 
decrease of conductance with the wire length for even C- 
wires in the anti parallel case is due to the presence of 
cr-bonds in these systems, which can act as tunnel bar- 
riers for electron conduction. In fact, recent calculations 
on (T-bonded structures have shown that the tunnel bar- 
rier increases with the increase of wire length^^, leading 
to an exponential decay in the electronic conduction. In 
contrast, in odd C-atom wires, the 7r-orbitals are highly 
delocalized, providing pathways for electron transfer and 
consequently leading to higher conductance as seen in 
Fig. 4. As discussed above, the super- exchange effect 
vanishes for a wire containing 14 C-atoms. 
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FIG. 4. Conductance (in units of e^/h) evaluated at the 
Fermi energy for both parallel (ferro) and anti parallel (anti- 
ferro) spin configurations in the carbon wires, as a function 
of the number of carbon atoms in the wires, n. 
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To understand the oscillatory pattern in conductance, 
we have also calculated the MuUiken charges and spin 
densities at individual atoms for different wire lengths. 
The spin density at the Co atoms, i.e. the difference be- 
tween the number of spin up and spin down electrons, are 
shown in Fig. 3. We see that for both parallel and anti- 
parallel magnetization states, the spin density at the Co- 
atoms oscillates with the number of C-atoms in the wire. 
Also, the (7 — TT conjugated wires show a higher atomic 
spin density at Co than the 7r-conjugated wires. This is 
not surprising since, as noted above, the 7r-conjugated 
systems have a stronger delocalization of spin compared 
to that in the a — -k conjugated wires. We have also 
calculated the magneto conductance (AfC) according to 
Eq.(4). Fig. 5 summarizes the MC values as a function 

120 I 1 




2 4 6 8 10 12 14 16 
Number of carbon atoms in wire 

FIG. 5. Dependence of magneto conductance, MC, eq.(4), 
on the number of carbon atoms in the atomic wire, n. 

of number of C-atoms in the wires, with MC displayed 
as a percentage. We find an oscillatory behavior in the 
magneto conductance, with a maximum value of 100% 
change in magneto resistance between the parallel and 
anti parallel magnetization states for wires containing 12 
and 14 C-atoms. This huge change in resistance between 
the two magnetization states suggests potential useful 
applications of these nanoscale materials for molecular 
magneto-electronics. 



bers of C atoms and acetylenic for even numbers of C 
atoms. The spin-polarized conductance was calculated as 
a function of number of C-atoms in the wire in the zero 
bias limit using Landauer's formalism. These length de- 
pendent calculations reveal an oscillatory pattern in con- 
ductance, with a significantly higher conductance arising 
for the parallel magnetization state compared to that for 
the anti-parallel magnetization state. The ground state 
C-wire structures containing up to 13 carbon atoms are 
found to have an anti-ferro magnetic spin configuration 
of the terminal Co atoms. In contrast, the wire with 14 
carbon atoms is found to have a parallel Co spin config- 
uration in the ground state. The energy difference be- 
tween the parallel and anti parallel magnetization states 
is found to be larger than fc^T at room temperature, sug- 
gesting that these two magnetization states are not in- 
terchangeable at normal operational temperatures. The 
stability of the anti ferromagnetic spin configuration be- 
tween the terminal Co atoms is seen to arise because of 
a super-exchange interaction that is facilitated by strong 
orbital overlap of the terminal magnetic Co atoms and 
the non-magnetic C-atoms of the wire. This effect van- 
ishes for the C-wire containing 14 C-atoms explaining the 
switch to a more stable parallel spin configuration. For 
carbon wires containing 12 and 14 C-atoms, we found 
almost no conductance in anti-parallel spin configura- 
tions, suggesting that the characteristic length for super- 
exchange interaction in these cr — tt conjugated carbon 
wires is about 20 A. A maximum value of 100% change 
in magneto resistance was obtained for carbon wires con- 
taining 12 and 14 carbon atoms. 
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IV. CONCLUSIONS 

We have investigated the chain length dependent mag- 
netic structures and energetics associated with highly 
conjugated C- wires sandwiched between magnetic Co 
atom species using the gradient corrected density func- 
tional approach. The Co-terminated wires show an al- 
ternation of structure between cumulenic for odd num- 
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